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Abstract 

Measurements of the production of jets of particles in association with a Z boson in pp collisions at 
y/s = 7 TeV are presented, using data corresponding to an integrated luminosity of 4.6 fb~^ collected 
by the ATLAS experiment at the Large Hadron Collider. Inclusive and differential jet cross sections in Z 
events, with Z decaying into electron or muon pairs, are measured for jets with transverse momentum 
PT > 30 GeV and rapidity \y\ < 4.4. The results are compared to next-to-leading-order perturbative 
QCD calculations, and to predictions from different Monte Carlo generators based on leading-order 
and next-to-leading-order matrix elements supplemented by parton showers. 



Prepared for submission to JHEP 



Measurement of the production cross section of 
jets in association with a Z boson in pp collisions 
at v^ = TTeV with the ATLAS detector 



Abstract: Measurements of the production of jets of particles in association with a Z bo- 
son in pp colUsions at \fs = 7 TeV are presented, using data corresponding to an integrated 
luminosity of 4.6 fb~^ collected by the ATLAS experiment at the Large Hadron Collider. 
Inclusive and differential jet cross sections in Z events, with Z decaying into electron or 
muon pairs, are measured for jets with transverse momentum px > 30 GeV and rapidity 
\y\ < 4.4. The results are compared to next-to-leading-order perturbative QCD calcula- 
tions, and to predictions from different Monte Carlo generators based on leading-order and 
next-to-leading-order matrix elements supplemented by parton showers. 



Contents 

1 Introduction 1 

2 Experimental setup 2 

3 Monte Carlo simulation 4 

4 Event selection 5 

5 Background estimation 6 

6 Detector-level results 7 

7 Correction for detector effects and combination of channels 8 

8 Systematic uncertainties 11 

9 Theoretical predictions 12 

10 Results and discussion 13 

10.1 Jet multiplicities 14 

10.2 Jet transverse momentum 16 

10.3 Angular distributions 19 

10.4 Distributions after VBF preselection 20 

10.5 Inclusive quantities 21 

11 Conclusions 23 

12 Acknowledgements 26 



1 Introduction 

The production of jets of particles in association with a Z boson^ at hadron colliders 
provides an important test of perturbative quantum chromodynamics (pQCD). Such events 
also constitute a non-negligible background for studies of the Higgs boson candidate [1, 2] 
and searches for new phenomena. In these searches, the multiplicity and kinematics of 
jets in Z + jets events are exploited to achieve a separation of signal from background. 
This procedure often introduces scales larger than the mass of the Z boson, resulting in 
large logarithmic contributions in the calculation of higher-order QCD corrections to the 
predicted Z + jets cross section [3, 4]. The measured Z + jets cross section can be compared 



^The notation Z refers to the complete Zj^* interference. 



directly to fixed-order predictions at next-to-leading-order (NLO) in pQCD [5-7] and to 
Monte Carlo (MC) generators based on next-to-leading-order or leading-order (LO) matrix 
elements supplemented by parton showers [8-10]. The simulations based on LO matrix 
elements are affected by large uncertainties in the factorization and renormalization scales 
and need to be tuned and validated using data. 

Measurements of the Z+jets cross section have been reported for lower jet energies 
and lower jet multiplicities in proton-antiproton collisions at a center-of-mass energy of 
^/s = 1.96 TeV [11-13] and in proton-proton collisions based on a data set of 0.036 fb~^ 
collected at ^/s = 7 TeV [14, 15]. This article extends these measurements, using 4.6 fb~^ 
of proton-proton collision data collected by the ATLAS experiment in 2011 at ^/s = 7 TeV. 
The large data set allows cross sections to be measured for the production of up to seven 
jets in association with a Z boson. Differential jet cross sections are accessible for large jet 
multiplicities and for energy regimes up to 1 TeV, which allows the modelling of the Z+jets 
process to be probed for typical phase-space regimes expected from new phenomena and 
from Higgs boson production, for example via vector-boson-fusion (VBF). 

Selected events contain a Z boson decaying into a pair of electrons or muons. As- 
sociated jets are identified in a rapidity (y-'"^*) range of \y^'^^\ < 4.4 and with transverse 
momentum {p!^ ) oi p!^ > 30 GeV. The measurements comprise inclusive and exclusive 
jet multiplicities for different phase-space constraints and differential jet cross sections as 
a function of the transverse momentum and the rapidity of the four jets with the largest 
transverse momentum ('leading jets'). Cross sections for events with at least two jets in 
the final state are measured as a function of the invariant mass {m^^) and the angular 
separation of the two leading jets. Differential cross sections in events with at least one jet 
are measured as a function of the scalar j?t sum of the jets (S't), of the scalar px sum of 
the leptons and jets {Ht), and the transverse momentum of the Z boson candidate {p^). 
The results of the measurements are unfolded for detector effects and quoted at the par- 
ticle (hadron) level, where they are compared to predictions from fixed-order NLO pQCD 
programs and from several MC generators. 

The paper is organized as follows. The detector and the data sample are described in 
the next section. Section 3 provides details of the simulations used in the measurements, 
while section 4 describes the lepton and jet reconstruction and the event selection. The 
estimation of background contributions is described in section 5 and selected uncorrected 
distributions are presented in section 6. The procedures used to unfold the measurements 
for detector effects and to combine electron and muon channels are detailed in section 
7. Systematic uncertainties are discussed in section 8. The NLO pQCD predictions are 
described in section 9. Measured cross sections are presented in section 10 and compared 
to generator and NLO pQCD predictions. Finally, section 11 provides a summary. 

2 Experimental setup 

The ATLAS detector [16] at the LHC covers nearly the entire solid angle around the colli- 
sion point. It consists of an inner tracking detector surrounded by a thin superconducting 
solenoid, followed by electromagnetic and hadronic calorimeters and a muon spectrometer 



incorporating three large superconducting toroid magnets (each with eight coils). The in- 
ner detector (ID) is immersed in a 2 T axial magnetic field and provides charged-particle 
tracking in the pseudorapidity^ range |r/| < 2.5. The high-granularity silicon pixel detector 
covers the vertex region and typically provides three measurements per track, the first hit 
being normally in the innermost layer. It is followed by the silicon microstrip tracker, which 
provides typically eight measurements (four space-points) per track. These silicon detectors 
are complemented by the transition radiation tracker, which covers a region up to \7]\ = 2.0. 
The transition radiation tracker also provides electron identification information based on 
the fraction of hits above a high energy-deposit threshold corresponding to transition radi- 
ation. The calorimeter system covers the pseudorapidity range \r]\ < 4.9. Within the region 
\r]\ < 3.2, electromagnetic calorimetry is provided by barrel and endcap high-granularity 
lead/liquid-argon (LAr) calorimeters. An additional thin LAr presampler covers |t/| < 1.8 
to correct for energy loss in material upstream of the calorimeters. Hadronic calorimetry is 
provided by a steel/scintillating-tile calorimeter, segmented radially into three barrel struc- 
tures within \r]\ < 1.7, and two copper/LAr hadronic endcap calorimeters, that cover the 
region 1.5 < \r]\ < 3.2. The solid angle coverage is completed in the region of 3.1 < |?7| < 4.9 
with forward copper/LAr and tungsten/LAr calorimeter modules optimized for electromag- 
netic and hadronic measurements respectively. The muon spectrometer (MS) comprises 
separate trigger and high-precision tracking chambers measuring the deflection of muons 
in a magnetic field generated by superconducting air-core toroids. The precision chamber 
system covers the region |7;| < 2.7 with three layers of monitored drift tubes, complemented 
by cathode strip chambers in the forward region, where the background is highest. The 
muon trigger system covers the range \7]\ < 2.4 with resistive plate chambers in the barrel, 
and thin gap chambers in the endcap regions. A three-level trigger system is used to select 
interesting events. The Level- 1 trigger is implemented in hardware and uses a subset of 
detector information to reduce the event rate to a design value of at most 75 kHz. This is 
followed by two software-based trigger levels which together reduce the event rate to about 
400 Hz. 

The analysis is based on a sample of proton-proton collisions at y/s = 7 TeV, collected 
in 2011 during periods of stable beam operation. Di-electron final states are selected with 
a trigger requiring at least two electrons of px > 12 GeV, using an electron identification 
similar to the one used in offline selection. Di-muon final states are selected with a trigger 
requiring at least one muon of pT > 18 GeV, using a higher-level trigger algorithm similar 
to the one used in the offline selection. The integrated luminosity used in both channels is 
4.64±0.08 fb"^ [17]. 



^ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The a;-axis points from the IP to the centre 
of the LHC ring, and the j/-axis points upward. Cylindrical coordinates (r, <j!>) are used in the transverse 
plane, <j) being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 8 as rj — ~ lntan(^/2). 



3 Monte Carlo simulation 



Monte Carlo event samples are used to determine background contributions, correct the 
measurements for detector effects, correct the theory calculations for non-perturbative ef- 
fects, calculate acceptance corrections, and estimate systematic uncertainties on the final 
results. 

Signal events (^ (— )• /U/x)+jets and Z (^ ee)+jets) are generated using ALPGEN 
v2.13 [8] interfaced to HERWIG v6.520 [18] for parton shower and fragmentation and to 
JIMMY v4.31 [19] for modelling interactions of the proton remnants, referred to as 'un- 
derlying event' in the following, using the AUET2-CTEQ61L tune [20]. In the following 
sections, the expression 'ALPGEN' refers to this version unless stated otherwise. Similar 
samples are produced with ALPGEN v2.14 interfaced to PYTHIA v6.425 [21] using the PE- 
RUGIA2011C [22] tune. For both ALPGEN samples, CTEQ61L [23] parton distribution 
functions (PDFs) are employed. Signal samples are also generated with SHERPA vl.4.1 
using the MEnloPS approach [10] and with MC@NLO v4.01 [24], interfaced to HERWIG, 
both using the CTIO [25] PDF set. The program PHOTOS [26] is used to simulate QED 
final state radiation (FSR) in the ALPGEN samples. QED-FSR simulation in SHERPA is 
based on the YFS method [27]. ALPGEN and SHERPA matrix elements are generated for 
up to five partons. The signal samples do not include Z+jets events produced via VBF. 
Based on generator-level studies, the expected contribution of these events to the measured 
cross sections is at the per-mille to per-cent level for the selections and kinematic ranges 
explored in this paper and always significantly below the statistical and systematic precision 
of the measurement. 

Background samples from M^ + jets and Z (— )• tt) +jets final states are generated sim- 
ilarly to the signal samples, using ALPGEN interfaced to HERWIG. The Vl^ + jets and 
Z + jets samples are normalized globally to next-to-next-to-leading-order (NNLO) pQCD in- 
clusive Drell-Yan predictions as determined by the FEWZ [28] program using the 
MSTW2008NNLO PDF set [29]. The uncertainties of about 5% are taken from an en- 
velope of predictions using different PDF sets and factorization and renormalization scales, 
as described in ref. [30]. Single-top-quark events are produced with AcerMC [31], inter- 
faced to PYTHIA, using CTEQ61L PDFs. Diboson processes {WW, WZ and ZZ) are 
simulated with HERWIG using the AUET2-L0* tune [20]. Reference cross sections for 
single-top-quark and diboson processes are calculated using the MC@NLO generator with 
the MSTW2008 PDF set [29]. The ti samples used for the relative normalization of final 
states in top-quark pair-production are generated with MC@NLO interfaced to HERWIG 
and with POWHEG [32, 33] interfaced to PYTHIA, both using the CTIO PDF set. 

All samples are processed through the GEANT4-based simulation [34, 35] of the AT- 
LAS detector. The simulation includes the modelling of additional pp interactions in the 
same and neighbouring bunch crossings (pile-up), with an average of nine interactions per 
crossing, that matches the distribution of interactions per crossing measured in data. 



4 Event selection 



Table 1 summarizes the kinematic regions in which Z bosons and jets are selected. They are 
defined to provide a good experimental coverage for the reconstruction of electrons, muons 
and jets in the event. Events with less than three tracks associated to the hard scattering 
vertex, defined as the vertex with the highest pT sum of its associated tracks, are discarded. 

Electrons are reconstructed from clusters of energy in the electromagnetic calorimeter 
matched to inner detector tracks. The electron candidates must have pT > 20 GeV and 
\r]\ < 2.47, excluding the transition region 1.37 < \rj\ < 1.52 between barrel and endcap 
electromagnetic calorimeter sections, and pass the 'medium' identification criteria described 
in ref. [36], re-optimized for 2011 conditions. No additional isolation requirement is applied, 
since non-isolated electron candidates are already suppressed by the identification criteria. 
Muon candidates are identified as tracks in the inner detector matched and combined with 
track segments in the muon spectrometer [37]. They are required to have pT > 20 GeV 
and \r]\ < 2.4. In order to achieve a sufficient rejection of multi-jet events, muons are re- 
quired to be isolated: the scalar sum of the transverse momenta of tracks within a cone 
of AR = y/{A(j))'^ -\- (Ar/)^ = 0.2 around the muon candidate must be less than 10% of 
the transverse momentum of the muon. All lepton pairs are required to have a separa- 
tion of AR > 0.2. The Z candidates are selected by requiring exactly two oppositely 
charged leptons of the same flavour. Their invariant mass [m ) must be within the range 
66 GeV < m^^ < 116 GeV. With this selection, 1228767 Z {-^ ee) and 1678500 Z {-^ ^u/i) 
candidate events are identified. 

Jets are reconstructed using the anti-kt algorithm [38] with a distance parameter 
R = 0.4. The inputs to the jet algorithm are topological clusters of energy in the calorime- 
ter [39]. The energies and directions of reconstructed jets in data and simulated events are 
corrected for the presence of additional proton-proton interactions, the position of the pri- 
mary interaction vertex, the measurement biases induced by calorimeter non-compensation, 
additional dead material, and out-of-cone effects, using detector simulation and a combi- 
nation of in-situ methods [39, 40]. Jets are required to have a transverse momentum above 
30 GeV and a rapidity of \y^^^\ < 4.4. Jets closer than 0.5 in AR to a selected lepton are 
removed. In order to reject jets from additional proton-proton interactions, the 'jet vertex 
fraction' is used. This is defined as the px sum of the tracks associated to the jet which 
are consistent with originating from the primary vertex divided by the pT sum of all tracks 
associated to the jet. The jet vertex fraction is required to be greater than 0.75 for jets with 
\r]\ < 2.4. The residual impact of additional proton-proton interactions on the distribution 
of the jet observables has been checked to be correctly simulated such that the unfolded 
cross sections are expected to be independent of the number of additional interactions. 
With this definition, 191566 Z (— )• ee) and 257169 Z (— t- ^fi) candidate events are selected 
with at least one jet in the final state. 





Z(^ee) 


^Hw) 


lepton pt 
lepton \r]\ 


PT > 20 GeV 

I77I < 1.37 or 1.52 < \r]\ < 2.47 


PT > 20 GeV 

|r/| < 2.4 


lepton charges 

lepton separation AR^^ 

lepton invariant mass m^^ 


opposite charge 

Ai?« > 0.2 

66 GeV < m^^ < 116 GeV 


jet pT 

jet rapidity yj°* 

lepton-jet separation AR^ 


p!^^ > 30 GeV 

|yjet| < 4_4 

AR^^ > 0.5 



Table 1. Summary of Z (— > ££) and jet selection criteria. 



5 Background estimation 

The selected data sample is expected to contain background events with two isolated leptons 
(it, diboson and Z {-^ tt) events), with one isolated lepton (W — )• ev, W ^ fiv and 
single-top-quark production) and without isolated leptons (multi-jet events). The total 
expected background fraction increases with the jet multiplicity (A'jet) from 2% (A'^jct > 1) 
to 20% (iVjet > 6). It is dominated by multi-jet processes, tt and diboson events for 
Z (+ < 1 jet) and by tt for larger jet multiplicities. The background is estimated using 
simulated samples, with the exception of the multi-jet and tt background contributions, 
which are derived from data. For these data-driven background estimates, the shape of the 
background contribution to each of the measured distributions is derived from a dedicated 
background-enriched sample in data. The background-enriched samples have been selected 
and normalized as described below. 

The multi-jet background contribution in the Z (— )• ee) +jets channel is estimated using 
a multi-jet enriched data template with two electron candidates which both pass a 'loose' 
selection but fail to pass the medium identification requirements [36] . The dedicated trigger 
used for the selection of this sample requires two clusters of energy in the electromagnetic 
calorimeter with pT > 20 GeV. This sample is dominated by jets misidentified as electrons 
in the final state. The normalization of this sample to the multi-jet background expected 
with medium requirements is extracted from a template fit in the invariant mass distribution 
for medium electrons {rrf^) as follows: A single combined fit is performed of the multi-jet 
template and the standard simulated signal and non-multi-jet background templates to the 
measured spectrum of the invariant mass for medium electrons in the extended mass range 
50 GeV < rrf^ < 150 GeV in the inclusive selection. Systematic uncertainties are assessed 
by varying the mass range and the binning in the fit, by using a different generator (SHERPA 
instead of ALPGEN) for the signal template, by varying the electron energy scale and 
resolution in the simulation and by allowing for a modification of the shape of the mass 
distribution in the multi-jet enriched sample. The multi-jet background to the measured 
inclusive jet multiplicities varies between (0.65 it 0.23)% for A'jet ^ 1 and (1.20 it 0.44)% 
for Ajet > 6. 

In the Z (— )• /i/i) + jets channel, heavy flavour production (with muons originating from 



b- and c-quark decays) and decays-in-flight of pions and kaons are tlie primary source of 
the multi-jet background, which is highly suppressed by the isolation requirement applied 
to the muon candidates. The multi-jet template is derived from a data sample where both 
muons fail the isolation requirement. The normalization factor is obtained by fitting the 
multi-jet template together with a template composed of the simulated signal and the non- 
multijet background events that pass the signal selection to the spectrum of the invariant 
mass of isolated muons (m'^^) measured in data in the range 40 GeV < m,'^^ < 150 GeV. In 
contrast to the Z (— )■ ee) + jets channel, the creation of the template and the normalization 
is performed separately for A^jet > 0, N^et > 1 and A'jet > 2. The normalization factor 
derived for A'jet > 2 is used for all higher jet multiplicities. The systematic uncertainty is 
assessed by replacing the multi-jet template with one formed from muons passing a loose 
isolation cut but failing the tight cut used to select signal muons. Multi-jet fractions vary 
between (0.25 ± 0.04)% for iVjet > 1 and (2.2 ± 2.2)% for iVjet > 6. 

The it background contributions in the Z (— )■ U) + jets samples are dominated by events 
where both W bosons decay leptonically. Since the kinematic properties of the jets in the 
final state are independent of the flavours of the two leptons, final states with one electron 
and one muon can be used to model the tt background contributions to Z (— )• ee) and 
Z (— )• fin) selections. The M-enriched sample is selected from data in the e^/i^ final state 
with kinematic requirements analogous to the Z (-^ ii)+ jets selection. The dedicated 
trigger used for the selection of this sample requires an electron with p^ > 10 GeV and 
a muon with pT > 6 GeV. For each of the observables, the number of W^+jets, Z+jets 
and diboson events expected from simulation in the tt-enriched sample is subtracted. The 
normalization from the e /i^ to the e^e~ and ^J.~^fi~ final states is calculated from ti samples 
generated with MC@NLO+HERWIG and with POWHEG+PYTHIA, separately for each 
jet multiplicity. Systematic uncertainties on the normalization arise from the choice of the 
generator, uncertainty on the lepton trigger, reconstruction and identification efficiency (see 
section 8) and on the electroweak background subtraction. The ti fractions vary between 
(0.80±0.05)% for Z (^ ee) + > 1 jet and (18.6±7.0)% for Z (^ ee) + > 6 jets and between 
(0.74 ± 0.03)% for Z {^ nfi) + > I jet and (18.1 ± 5.3)% for Z (^ ^/i) + > 6 jets. 

6 Detector-level results 

Measured and expected distributions of the jet observables have been compared at the re- 
construction level, separately in the electron and muon channels. As an example, figure 1 
shows the dilepton invariant mass in events with at least one jet in the final state, as well as 
the inclusive jet multiplicity. For the signal, both ALPGEN and SHERPA expectations are 
shown. In this figure, W — t- eu, Z (— >• tt) and diboson processes are summarized as 'elec- 
troweak' background and tt and single-top processes are referred to as 'top' background. 
For figures 1(a) and 1(b), the selection has exceptionally been extended beyond the fidu- 
cial invariant mass range, in order to demonstrate in addition the reasonable agreement 
between data and expectations for dilepton mass sideband regimes with larger background 
fractions. Table 2 shows, for the electron and muon channels separately, the observed num- 
ber of events for the different jet multiplicities in the final state compared to expectations 
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Table 2. Numbers of events expected and observed in data that pass the Z {—>■ ee)+jets and 
Z (— >■ nfj.) + jets selections as a function of the inclusive jet multiplicity. The expected numbers are 
rounded according to the combined statistical and systematic uncertainty. ALPGEN has been used 
to simulate the signal events. 

for signal (ALPGEN) and background processes. The combined statistical and systematic 
uncertainties on the total expectation increases from 6% to 30% with increasing jet multi- 
pHcity. The data are consistent with predictions by the generators ALPGEN and SHERPA, 
which gives confidence that the simulated samples, which are used in the unfolding, provide 
a reasonable description of the event kinematics and of the detector response. 

7 Correction for detector effects and combination of channels 



The cross sections in this article are quoted at the particle level, which corresponds to 
'dressed' muons and electrons, calculated using final-state leptons from the Z decay for 
which collinear radiation in a cone of AR < 0.1 is added to the lepton four-momentum. 
Particle jets are clustered from all final-state particles (decay length cr > 10 mm) excluding 
the dressed Z decay products. The phase-space requirements are the same as in the selection 
at reconstruction level (see table 1). 

After subtracting the expected background contributions, the data distributions in 
each channel are unfolded to the particle level using an iterative technique [41]. Response 
matrices are calculated for each observable, using Z + jets samples generated with ALPGEN. 
Before entering the iterative process, the data are corrected for the fraction of reconstructed 
events in the ALPGEN sample which do not match to a particle-level equivalent. The 
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Figure 1. Numbers of events observed in data and predicted in simulation that pass the 
Z (— ^ ee)+jets and .^ (— >■ ////)+ jets selection as a function of the invariant mass of the Z can- 
didate, (a) vff^ and (b) m'^'^, for events with at least one jet with pi^ > 30 GeV and |yj®*| < 4.4, 
and as a function of the inclusive jet multiplicity, -A'jet, in (c) di-electron and (d) di-muon events. 
The individual contributions of the various backgrounds are also shown, as detailed in the legend. 
The hatched band corresponds to the combined statistical and systematic uncertainty on the pre- 
diction, obtained using ALPGEN to model the Z + jets process. The error bars on each data point 
show the statistical uncertainty. The bottom panel shows the corresponding MC/data ratio. The 
shaded band corresponds to the total systematic uncertainty and the error bars to the statistical 
uncertainty on the MC/data ratio. 



number of iterations, typically two or three, is optimized for each observable using a ~)(^ 
comparison of generated and unfolded reconstructed Z + jets events from the generators 



SHERPA and MC@NLO. 

The uncertainties from the hmited number of events in data are propagated into 
the particle-level cross sections using a Monte Carlo method. One thousand pseudo- 
experimental spectra are generated by fluctuating the content of each bin according to 
the statistical uncertainty. The unfolding procedure is applied to each pseudo-experiment, 
and the r.m.s. of the results is taken as the statistical uncertainty. Systematic uncertainties 
arising from the unfolding procedure are estimated by comparing with an iterative un- 
folding based on response matrices and corrections derived from SHERPA. The statistical 
uncertainties of the response matrices are propagated into systematic uncertainties on the 
unfolded cross sections using pseudo-experiments. 

The cross sections measured in the electron and muon channels are extrapolated to a 
common phase-space region, derived from table 1 by extending the rj range of the leptons to 
l^iepj ^ 2.5^ using global acceptance corrections derived from ALPGEN Z + jets Monte Carlo 
samples, reweighted to the CTIO PDF set. The corrections are of the order of 14% and 5% 
for the electron and muon channel, respectively. Systematic uncertainties are estimated by 
comparing with corrections obtained using the corresponding SHERPA Z+jets sample and 
the original ALPGEN sample. Total uncertainties on the corrections are calculated as the 
quadratic sum of the statistical and systematic uncertainties and amount to 0.2-0.3%. The 
extrapolated cross sections measured in the electron and muon channels are in agreement. 

For each observable, the extrapolated cross sections are combined using the averaging 
procedure introduced in ref. [42], which accounts for systematic uncertainties (bin-to-bin 
correlated and uncorrelated) proportional to the central values of the respective cross sec- 
tions. The weights of the individual cross-section measurements (;U^) in channel k (ee or 
/Xyu) and bin i in the combined cross sections (?7i* ) are derived by minimizing the following 
X^ function [42]: 

where bj denote the shift introduced by a correlated systematic error source j normalized 
to its respective standard deviation. The relative statistical and uncorrelated systematic 
uncertainties on fi^, are denoted by 61^^^ j^ and 61^^^^^. ^ and the variable 7* j^ quantifies the 
influence of the correlated systematic error source j on the measurement i in the channel 
k. 

The following bin-to-bin correlated systematic sources are taken into account: nor- 
malization of the multi-jet background, lepton energy scale and resolution, lepton recon- 
struction, identification and trigger efficiencies and normalization of tt, electroweak and 
single-top background contributions, the latter three treated as correlated between the 
channels. Bin-to-bin correlated systematic sources which have the same impact in both 
channels do not enter in the combination procedure. These are the individual components 
of the jet energy scale, the jet energy resolution, the luminosity, the unfolding procedure, 
and the extrapolation factor. The uncertainties from these sources on the combined result 
are taken as the weighted average of the corresponding uncertainties on the electron and 
muon measurements. 
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8 Systematic uncertainties 

The kinematic ranges and the binning are chosen such that the statistical uncertainty of 
the measurement is comparable to or smaller than the systematic uncertainty. The relative 
systematic uncertainties on the cross sections measured in each channel are derived for 
each observable by propagating systematic shifts from a set of independent sources through 
the response matrices and the subtracted background contributions into the unfolded data. 
The resulting systematic uncertainties for each source in each channel are symmetrized in 
order to mitigate the impact of statistical fluctuations and are combined in the averaging 
procedure. 

The uncertainty on the jet energy scale (JES), determined from the combination of 
methods based on MC and in-situ techniques used to determine the scale, constitutes the 
dominant component of the total systematic uncertainty. It is propagated through the 
analysis using 14 independent components fully correlated in pj, [39, 40]. They account for 
uncertainties on the different in-situ measurements which enter the jet calibration, on the 
jet flavour and on the impact of pile-up and close-by jets. The uncertainty on the jet energy 
resolution, derived from a comparison of the resolution obtained in data and in simulated 
dijet events, is propagated into the final cross section by varying the energy resolution of 
the simulated jets. Uncertainties on the normalization of the background expectations, for 
simulated and data-driven background contributions respectively, are treated as correlated 
between bins and are propagated to the measured cross sections by unfolding the data 
distributions after the subtraction of the systematically shifted background. The statistical 
uncertainties of the background contributions are added quadratically to the statistical 
uncertainties of the data. The uncertainty from the unfolding process is derived from 
the different components discussed in section 7, which are considered to be uncorrelated. 
Systematic uncertainties on electron and muon trigger efficiencies, energy scale, resolution, 
reconstruction and identification efficiencies are derived from the comparison of tag-and- 
probe results in data and simulated events [36, 37]. 

Table 3 summarizes the systematic uncertainties on the Z + jets cross sections as a func- 
tion of the inclusive jet multiplicity and of p!^ of the leading jet separately for the electron 
and muon channels. The uncertainty on the integrated luminosity of 1.8% translates into 
comparable uncertainties on the measured cross sections. The total uncertainties on the 
inclusive jet cross sections range from 8% for A^jgt > 1 to 16-17% for A'jet > 4, dominated 
by the JES uncertainty. 

The uncertainty on cross-section ratios, -R>(„+i)/>n'^ ) for successive jet multiplicities n 
is significantly reduced due to the strong correlations between the lepton and jet reconstruc- 
tion and calibration uncertainties in neighbouring jet bins and amounts to a total of 3-4% 
for R>2/>1 ^i^d higher multiplicities, which are of interest in this article, dominated by the 
residual JES uncertainty. The large JES uncertainties in the forward region propagate into 
uncertainties on the unfolded cross sections at the level of 20% (30%) for jet rapidities of 
lyJ^^t] = 3.0 (4.0). This is reflected in large jet energy scale uncertainties on the cross section 



^For simplicity, n is used in the subscript instead of A^jot 
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Z (^ ee) 


>ljet 


> 2 jets 


> 3 jets 


> 4 jets 


^^ in [30-500 GeV] 


electron reconstruction 
jet energy scale, resol. 
backgrounds 
unfolding 


2.8% 

7.4% 

0.26% 

0.22% 


2.8% 
10.1% 
0.34% 
0.94% 


2.8% 

13% 

0.44% 

1.2% 


2.8% 

17% 

0.50% 

1.9% 


2.6-2.9% 
4.3-9.0% 
0.2-3.2% 
1.4-6.8% 


total 


7.9% 


10.5% 


13% 


17% 


5.5-12.0% 


Z (^ ii[i) 


>ljet 


> 2 jets 


> 3 jets 


> 4 jets 


^^ in [30-500 GeV] 


muon reconstruction 
jet energy scale, resol. 
backgrounds 
unfolding 


0.86% 

7.5% 

0.093% 

0.30% 


0.87% 
9.9% 
0.20% 
0.68% 


0.87% 

13% 

0.41% 

0.52% 


0.88% 
16% 

0.66% 
1.3% 


0.8-1.0% 
3.2-8.7% 
0.1-1.9% 
0.5-6.2% 


total 


7.6% 


10.0% 


13% 


16% 


4.4-10.2% 



Table 3. Systematic uncertainties on the cross sections for Z (— )■ ee) + jets and Z (— )■ \i\i) +jets as a 
function of the inclusive jet multiplicity and as a function of the transverse momentum, pi^ , of the 
leading jet for events with at least one jet with 'p'^ > 30 GeV and ly-'^*] < 4.4. The rows labelled 
'electron reconstruction' and 'muon reconstruction' include uncertainties on trigger, reconstruction 
and identification, energy scale and resolution. 

for events with large rapidity distance (lAy-'-'l) between the leading jets, which combine 
with the unfolding uncertainties to total uncertainties of 20% (50%) for jAy-'-'l = 3.0 (4.0). 



9 Theoretical predictions 

Fixed-order calculations at NLO pQCD for the production of Z (+ > 1 jet) up to 
Z {+>4: jets) are computed using the BlackHat+SHERPA program [5-7]. CTIO PDFs [25] 
are employed and renormalization and factorization scales are set to Ht^/2, where Ht^ is de- 
fined event-by-event as the scalar sum of the pT of all stable particles/partons. The anti-A;^ 
algorithm with R = 0.4 is used to reconstruct jets at the parton level. Systematic uncer- 
tainties on the predictions related to PDF uncertainties are computed from the 52 CTIO 
eigenvectors at 68% confidence level [25]. The uncertainties on the cross sections increase 
from 1% for (iVjet > 1) to 3% for (iVjet > 4) and from 1% to 5% with pj^* of the leading jet 
between 30 GeV and 500 GeV. Additional changes in the PDFs due to the variation of the 
input value for the strong coupling constant as at the Z-boson mass scale by ±0.001 around 
its nominal value as{mz) = 0.118 introduce uncertainties on the predicted cross sections in 
the range of 1% to 3% for Z (+ (1 — 4) jets). These are added in quadrature to the PDF 
uncertainties. Scale uncertainties are estimated by variations of the renormalization and 
factorization scales to one half and two times the nominal scale. The scale uncertainties for 
different parton multiplicities are assumed to be uncorrelated. For inclusive calculations, 
the scale variations translate into variations of the cross section by 4% to 13% as A'jet in- 
creases and by 2% to 18% with increasing p^^ of the leading jet. For exclusive final states, 
the scale uncertainties are calculated using the prescription of ref. [43]. For comparison, the 
theory /data ratios presented in section 10 also show the scale uncertainty resulting from a 
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simple variation of the renormalization and factorization scales by a factor of two, assuming 
the uncertainties to be correlated for different parton multiplicities. The scale uncertainties 
constitute the dominant uncertainties in most kinematic regions. 

The NLO fixed-order calculations at the parton level are corrected to the particle level 
for the underlying event and for effects of fragmentation and of QED final-state radiation 
(QED-FSR). Parton-to-hadron correction factors [6 ) approximately account for non- 
perturbative contributions from the underlying event and fragmentation into particles. For 
each observable, the correction factor is estimated using simulated Z + jets samples, pro- 
duced with ALPGEN with the HERWIG cluster fragmentation in which JIMMY models 
the underlying event using the AUET2-CTEQ61L [20] tune. It is calculated as the bin-by- 
bin ratio of the nominal distribution at the particle level to the one obtained by turning 
off both the interactions between proton remnants and the fragmentation in the simulated 
samples. The non-perturbative corrections are also computed using ALPGEN samples, 
this time interfaced to PYTHIA, where the correction corresponds to the combined effect 
of string fragmentation and of the underlying event predicted by the PERUGIA2011C [22] 
tune. The difference is taken as a systematic uncertainty. The combined nominal correction 
is 7% in the low p!^ region and decreases with increasing pl^ towards zero. The correction 
factors for the inclusive iVjet distributions are about 3-4%. Nonperturbative corrections for 
quantities calculated with several jets include implicitly the corrections for all jets. The 
statistical and the symmetrized systematic uncertainties on 5^'^'^ are added in quadrature 
to the total uncertainty from the BlackHat+SHERPA calculation. 

The QED-FSR correction factors (5^^^) are determined using Z+jets samples pro- 
duced with the ALPGEN generator, interfaced to PHOTOS [26], by calculating the ex- 
pected cross sections both with the lepton four-momentum before final-state photon radi- 
ation ('Born level'), and with dressed leptons. The correction factors are about 2% for the 
electron and muon channels. They do not show a significant A'jet dependence and are stable 
with respect to the jet rapidity and for large jet transverse momentum. Systematic uncer- 
tainties are derived by comparing with (5^^^ obtained using a Z + jets sample produced 
with the SHERPA generator [9] which generates QED-FSR using the YSF method [27]. 
The differences between the two predictions are usually at the per-mille level. 

10 Results and discussion 

For each observable, the spectrum measured in data is unfolded to the particle level. After 
extrapolation and combination of electron and muon channels, the results are compared 
with calculations from BlackHat+SHERPA, corrected to the particle level, and with 
predictions by ALPGEN, SHERPA and MC@NLO. Both ALPGEN and SHERPA em- 
ploy matrix elements for up to five partons. Higher multiplicities are generated by the 
parton shower. In contrast, MC@NLO generates the Drell-Yan process at NLO preci- 
sion, which includes the real emission of one additional parton. All higher parton multi- 
plicities are generated by the parton shower. Inclusive and differential cross sections for 
Z {-^ U)+ > njets are compared with BlackHat+SHERPA fixed-order pQCD calcula- 
tions for Z+ > n partons, which provide a NLO estimate for the respective parton multi- 
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plicity, including the real emission of one additional parton. Measured cross sections as a 
function of the jet multiplicity and their ratios are detailed in table 4. Tabulated values of 
all observed results are available in the Durham HEP database [44]. 

10.1 Jet multiplicities 

Figure 2(a) presents the absolute cross sections for inclusive jet multiplicities for up to 
seven hadronic jets in the final state. The ratios -R>(n+i)/>n of cross sections for two 
successive multiplicities, presented in figure 2(b), provide a more precise measurement of 
the QCD process, due to the cancellation of part of the systematic uncertainty. The data are 
consistent with BlackHat+SHERPA calculations and with predictions of the generators 
ALPGEN and SHERPA. The MC@NLO parton shower underestimates the observed rate 
for additional jet emission by a factor of two, which leads to large offsets to the data for 
higher jet multiplicities. For this reason, in subsequent figures the MC@NLO predictions 
are only shown for Z (^ £i)+ > 1 jet selections, where the parton corresponding to the 
NLO real emission can be expected to yield a reasonable description of the kinematics. 

Exclusive jet multiplicities at the LHC are expected to be described by means of two 
benchmark patterns, 'staircase scaling' with R(n+i)/n constant and 'Poisson scaling' with 
Rin+i)/n inversely proportional to n [3, 45], which provide limiting cases for certain kine- 
matic conditions. While for high multiplicities a flat exclusive jet multiplicity ratio is 
derived from the non-abelian nature of QCD FSR, at low multiplicity the jet multiplicity 
ratio is flat due to the combined effect of a Poisson-distributed multiplicity distribution and 
parton density suppression [3] . The emission of the first parton should be suppressed more 
strongly than the subsequent parton emissions. The underlying Poisson scaling is expected 
to emerge after introducing large scale differences between the core process (Z (+1 jet)) and 
the pl^ of the second leading jet. Two selections are chosen to test the two benchmark 
scenarios: (a) the standard Z+jets selection and (b) events where the leading jet has a 
transverse momentum in excess of 150 GeV. 

Figure 3(a) presents the ratios R(n+i)/n of cross sections for two successive exclusive 
multiplicities for the standard Z + jets selection. The comparitatively large scale uncertain- 
ties on the pQCD predictions result from the prescription of ref. [43], assuming the scale 
variations to be uncorrelated across the jet multiplicities. For comparison, the total un- 
certainty calculated using a naive scale variation, and a reduced uncertainty that does not 
include any scale uncertainty are also shown. The data are consistent with the central values 
of the BlackHat+SHERPA calculations and with predictions by the generators ALPGEN 
and SHERPA. The cross-section ratios show an approximately linear dependence on the jet 
multiplicity with a small slope. A linear fit -R(„+i) /„ = Rq + -^ • n oi the observed mul- 
tiplicity ratio starting with i^2/i yields Rq = 0.232 ± 0.009 and dR/dn = -0.011 ± 0.003. 
The uncertainties include a systematic contribution, derived from a series of fits to system- 
atic variations of the multiplicity ratio. The flat staircase pattern provides an acceptable 
approximation of the observed scaling behaviour for the standard Z+jets selection. The 
observation is consistent with results presented in [15] on the smaller data set collected in 
2010. 
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Figure 2. (a) Measured cross section for Z (— > U) +jets as a function of the inclusive jet multiplicity, 
iVjet, and (b) ratio of cross sections for successive inclusive jet multiplicities. The data are compared 
to NLO pQCD predictions from BlackHat+SHERPA corrected to the particle level, and the 
ALPGEN, SHERPA and MC@NLO event generators (see legend for details). The error bars indicate 
the statistical uncertainty on the data, and the hatched (shaded) bands the statistical and systematic 
uncertainties on data (prediction) added in quadrature. 

Figure 3(b) presents the exclusive jet multiplicity ratio for events where the leading 
jet has a transverse momentum in excess of 150 GeV. The observed ratio i2(„+iw„ is now 
steeply increasing towards low jet multiplicities, a pattern described by the central values of 
the BlackHat+SHERPA calculations, by the generator ALPGEN and approximately also 
by SHERPA. The observed cross-section ratios have been fitted with a pattern expected 
from a Poisson-distributed jet multiplicity with the expectation value n, i?(n+i)/n = ^- The 
Poisson scaling provides a good overall description of the jet multiplicity observed in data 
for the selected kinematic regime, with n = 1.02 ± 0.04, where the uncertainty includes 
statistical and systematic components. 

The scaling pattern is also investigated for a preselection typically employed in the 
selection of particles produced via vector boson fusion (VBF). Figure 4 presents the absolute 
cross section as a function of the exclusive jet multiplicity and R(n+i)/n after requiring two 
jets with m^^ > 350 GeV and jAy-'-'l > 3.0, in the following referred to as 'VBF preselection'. 
The data are consistent with the BlackHat+SHERPA prediction. SHERPA describes the 
multiplicity well whereas ALPGEN overestimates R^/2- 
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Figure 3. (a) Ratio of cross sections for successive exclusive jet multiplicities, A'jot, in events se- 
lected with the standard selection and (b) in events with at least one jet with •p'^ > 150 GeV and 
|yjet| < 4 4^ xhe data are compared to NLO pQCD predictions from BlackHat+SHERPA cor- 
rected to the particle level, and the ALPGEN, SHERPA and MC@NLO event generators (see legend 
for details) . The error bars indicate the statistical uncertainty on the data, and the hatched (shaded) 
bands the statistical and systematic uncertainties on data (prediction) added in quadrature. The 
shaded bands on the theory calculations show the systematic uncertainty excluding the scale un- 
certainty (dark shaded) and the total systematic uncertainties using the naive approach (medium 
shaded) and the nominal approach (light shaded) to derive the scale uncertainty (see section 9). 
The figures include (a) a linear fit R(n+i)/n = ^o -|- 4^ • n in the range i?2/i < ^(n-i-i)/n < ^5/4 
and (b) a Poisson fit i?(„_|_i)/„ = ^ to the data points, with the free parameters Ro, ^ and n. 



10.2 Jet transverse momentum 

Differential cross sections witli respect to tlie jet transverse momentum, p!^ , provide a test 
of pQCD over a large kinematic range. In particular, when p!^ exceeds the scale given by the 
gauge boson mass, NLO/LO K-factors can be large due to the presence of QCD corrections 
of the order of a^ln {pl^ /mz) [4]. In addition, higher-order electroweak corrections are 
expected to reduce the cross section with increasing transverse momentum of the Z boson 
candidate, by 5-20% for 100 GeV < p^4 < 500 GeV [46]. 

Figures 5 and 6 show the cross section as a function of p!^ of the first, the second, the 
third and the fourth leading jet (in descending order of pl^ ) for events with at least one, 
two, three and four jets in the final state, respectively. The cross sections are normalized to 
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Figure 4. (a) Measured cross section for Z {-^ £^)+jets as a function of the exclusive jet multi- 
plicity, -/Vjct, and (b) ratio of the cross sections for two successive multiplicities, in events passing 
the VBF preselection (at least two jets with jsi^* > 30 GeV and lyj''*] < 4.4 and m" > 350 GeV and 
|Aj/-'^| > 3.0 for the two leading jets). The other details are as in Figure 3. 



the inclusive Z (— >• ii) cross section, which reduces the systematic uncertainties connected 
to lepton identification and integrated luminosity. The fixed-order NLO predictions by 
BlackHat+SHERPA are consistent with the data for all jet multiplicities. 

For the leading jet, the precision of the measurement exceeds the precision of the theory 
prediction. While ALPGEN predictions for the pl^ spectrum of the second to fourth leading 
jet are consistent with the data, the p!^ spectrum of the leading jet is predicted to be too 
hard for larger values of p!^ . SHERPA is characterized by offsets to the data at the level of 
5-15%, consistent with the observations presented in figure 2(a) for the inclusive jet cross 
section. MC@NLO predicts a too soft pl^ spectrum, resulting in a discrepancy with the 
data by one order of magnitude for large p!^ . This is attributed to the fact that the fraction 
of events with a second resolved jet, which in MC@NLO is modelled via the parton shower, 
increases considerably with p!^ of the leading jet (see figures 3(a) and 3(b) for small and 
larger p!^ (leading jet)). A too soft p!^ spectrum of the parton shower will hence result in 
an increasing discrepancy between the MC@NLO prediction and the data. 

Figure 7(a) shows the cross section as a function of p^r^ of the leading jet, normalized 
to the inclusive Z (— )• ii) cross section, when a veto on a second jet is applied. A better 
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Figure 5. (a) Measured cross section for Z (— > U) +jets as a function of the transverse momentum, 
p^r^ , of the leading jet for events -with at least one jet -with -f}^ > 30 GeV and |j/J°*| < 4.4 in the 
final state and (b) as a function of p!^ of the second leading jet for events with at least t-wo jets. 



The cross sections are normalized to the inclusive Z (- 
in Figure 2. 



i£) cross section. The other details are as 



agreement between the predicted and observed cross-sections is observed. For events with 
at least two jets, figure 7(b) shows cross section as a function of the ^r^; ratio of the two 
leading jets, normalized to the inclusive Z (— >■ ii) cross section. ALPGEN overestimates the 
cross section for events with a pj, ratio of the leading jets in the range of 0.1-0.2. SHERPA 
underestimates the cross section as a function of the p!^ ratio by «15%, consistent with 
the results presented in figure 2(a). 

In a complementary approach, the cross section is measured as a function of the px 
of the recoiling Z boson, reconstructed from the momenta of the two leptons. The results 
are presented in figure 8 for both the inclusive and the exclusive Z (+1 jet) selection, nor- 
malized to the inclusive Z (— )■ ii) cross section. Both ALPGEN and SHERPA predict a 
too hard p^ spectrum, in particular in the inclusive case. The discrepancy with the data is 
comparable to the expected higher-order electroweak corrections [46] although higher-order 
QCD corrections could equally account for this. The BlackHat+SHERPA Z {+ > 1 jet) 
fixed-order calculation for the inclusive final state is too soft whereas for the exclusive final 
state the central predictions are closer to the observed spectrum. This result is attributed to 
missing higher jet multiplicities in the fixed-order calculation and will be discussed in more 
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Figure 6. (a) Measured cross section for Z (— > U) +jets as a function of the transverse momentum, 
j}^ , of the third leading jet for events with at least three jets with pl^ > 30 GeV and \y^'^^\ < 4.4 
in the final state and (b) as a function of p!^ of the fourth leading jet for events with at least four 



jets. The cross sections are normalized to the inclusive Z (- 
are as in Figure 2. 



££) cross section. The other details 



detail in section 10.5. The comparison with BlackHat+SHERPA yields no indication for 
missing higher-order electroweak corrections in the large-p^ region. Consistent with the 
results presented for the pj, spectrum of the leading jet, MC@NLO describes the exclusive 
Z (+1 jet) final state better than the corresponding inclusive final state. 

10.3 Angular distributions 

Figures 9 and 10 show the rapidity spectrum of the four leading jets, normalized to the 
inclusive Z (-^ U) cross section. Both BlackHat+SHERPA and SHERPA predict rapidity 
spectra for the leading jet that are somewhat wider than observed in the data. ALPGEN 
predictions are compatible with the measurements. 

Figure 11 presents the separation in rapidity, |A?/-'-'|, and the invariant mass, m^^ , of 
the two leading jets, normalized to the inclusive Z (— )• ii) cross section. The predictions by 
BlackHat+SHERPA and ALPGEN are consistent with the data. SHERPA overestimates 
the cross section for large lAy-'-'l, consistent with the too wide rapidity spectra. 

Differential jet cross sections as a function of angular distances {Acp^^ and AR^^) be- 
tween the two leading jets are presented in figures 12(a) and 12(b), respectively, normal- 
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Figure 7. (a) Measured cross section for Z (— > ££) + jets as a function of the jet transverse momen- 
tum, p!^ , for events with exactly one jet with pl^ > 30 GeV and |j/^°'| < 4.4 in the final state and 
(b) as a function of the ratio of p!^ of the second leading jet to p!^ of the leading jet for events 
with at least two jets. The cross sections are normalized to the inclusive Z (— )■ £i) cross section. 
The other details are as in Figure 3. 

ized to the inclusive Z (— )• ii) cross section. The azimuthal distance is well modelled by 
ALPGEN and by BlackHat+SHERPA. The tendencies observed in the modelling of the 
distance in (j) and in rapidity are reflected in the measurement of the AR spectrum of the 
leading jets. SHERPA models a too flat spectrum for both A(j) and AR. The offset of 
15% of the SHERPA prediction from the observed cross section in the bulk of the data in 
figures 11 and 12 is consistent with the results presented in figure 2(a) for the inclusive 
Z (—7- £i)+ >2 jets cross section. 

10.4 Distributions after VBF preselection 

A veto on a third jet is used to reject Z+jets background in selections of Higgs boson 
candidates produced by VBF. Figure 13 shows the transverse momentum and rapidity dis- 
tributions of the third jet after the VBF preselection, as defined in section 10.1, normalized 
to the inclusive Z (— )• U) cross section. The predictions by BlackHat+SHERPA, ALP- 
GEN and SHERPA are consistent with the measurements. Figure 14 shows the fraction of 
events which have fulfilled the requirements of a VBF preselection that pass in addition a 
veto on a third jet in the central region (|?7| < 2.4) as a function of the minimum trans- 



20 



10- 

10-= 

10-= 

10-' 

10-= 

10-' 

10' 
1.4 
1.2 

1 
0.1 
0.6 
1.4 
1.2 

1 
0.8 
0.6 
1.4 
1.2 

1 
0.1 
0.6 



ATLAS 

I L dt = 4.6 fb ' 

anti-k, jets, R = 0.4 

pij* > 30 GeV, l/'l < 4.4 



Z/y*(^ IT)+ > 1 jet (l=e,n) 

-;•- Data 201 1 (is = 7 leVn 

-m- ALPGEN 

-i- SHERPA 

^^ MC@NLO 

-T— BlackHat + SHERPA - 




- BlackHat + SHEBPA 



^--.~. .^.^i-^^mm 



i4ii^i" ' 



^ ^^^^^v/x//M^^'^ 



^ 



50 100 150 200 250 300 350 400 450 

p![ [GeV] 



^ 10-' 



^ 10- 

1 10-' 
10-= 
10-' 
10-' 



nj 



1.5 



0.5 



1.5 



0.5 



ATLAS 

L dt = 4.6 fb' 



anti-k, jets, R = 0.4 

<^^ > 30 GeV, l/'l < 4.4 



z/y'(^ rr)+i jet (i=e,n) 

,^»- Data 2011 (fi=7TeVr 

-■- ALPGEN 

-*- SHERPA 

^- MC@NLO 

-▼- BlackHat + SHERPA 



- BuickHat + SHERPA 



I I I I I I 



S -15- -"-ALPGEN 



0.5- 



- T ^ - f^ T - t --f - 



I I I I I I 



^^H-^ 



1 1 1 1 1 1 



1 '4'^4^^^ - ^*^^-^ ^ ^^-^^'^'iv'//y^-: 



50 100 150 200 250 300 350 400 450 

p; [GeV] 

(b) 



Figure 8. (a) Measured cross section for Z (— > ££) + jets as a function of the transverse momentum 
of the Z candidate, p^, in events -with at least one jet -with p!^ > 30 GeV and |yj°* | < 4.4 in the final 
state and (b) as a function of p^ in events with exactly one jet. The cross sections are normalized 
to the inclusive Z (— >■ £i) cross section. The other details are as in Figure 3. 

verse momentum of the veto jet, referred to as 'jet veto efficiency' in the folio-wing. The 
results are sho-wn at detector level, separately for the Z — )• ee and the Z — )■ fifi channel. 
The overestimate of Rs/2 in ALPGEN (see figure 4) leads to an underestimate of the veto 
efficiency, particularly for the \ow-pl^ regime. SHERPA predicts the veto efficiencies better. 



10.5 Inclusive quantities 

Quantities based on inclusive pT sums of final-state objects, such as Ht or St, are often 
employed in searches in order to enrich final states resulting from the decay of heavy parti- 
cles. Reference [47] reports a discrepancy between fixed-order pQCD calculations and data 
for moderate energy regimes in VF+jets events which can be mitigated by including higher 
jet multiplicities in the theoretical calculations by means of 'exclusive sums' [48]. 

Differential cross sections of Z (+ > 1 jet) events as a function of i?T and St, nor- 
malized to the inclusive Z (— )• U) cross section, are presented in figure 15. ALPGEN 
predicts slightly too hard spectra for both variables in line with the too hard spectrum 
for p}^ . SHERPA predictions show an offset of 10-15% to the data. The softer spectra 
from BlackHat+SHERPA, based on a Z (+ > 1 jet) fixed-order NLO calculation, deviate 



21 



ATLAS 

I L dt = 4.6 fb ' 

anti-k, jets, R = 0.4 

p"*' > 30 GeV, l/'l < 4.4 



Z/y*(^ IT)+ > 1 jet (l=e,n) - 
'-It- Data 2011 (is = 7TeV) 
-■- ALPGEN 
-i- SHERPA 

MC@NLO 

BlackHat + SHERPA " 




^ 10-'- 



ATLAS Z/Y'Hrr)+>2jets(l=e,n) 

I L dt = 4.6 fb' ^<j^ Data 2011 (S = 7 TeV) 

anti-k, jets, R = 0.4 -■- ALPGEN 

p^;' > 30 GeV, l/'l < 4.4 ^4- SHERPA 

— T- BlackHat + SHERPA 



^^*t. 



^^A 



"?*2 



lO'-' - 






1.4- 
1.2- 



i-H 



- BlackHat + SHERPA 



i-H 



1 wy^'y'^yW/Cr^y////^^/, 



0.8- 

^ ■'* -■- ALPGEN 

1.2h 

1 
0.8- 






| H ll| 



'^^'T^^'^^^/ f^ VVA ^/ M i i^^M^^^M 



0.6^ 
1.4^ 
1.2- 

0.8^*" 
0.6- , 



I I I I I I 
- SHERPA 

■/// • /////, 



| H ll| 



22:53?Ri35P^^ 



2 2.5 3 3.5 4 

ly'"*! (leading jet) 



0.5 



1.5 



2.5 



3 3.5 4 
ly'"*! (2nd leading jet) 



(b) 



Figure 9. (a) Measured cross section for Z (— > it) + jets as a function of the absolute value of the 
rapidity, |2/''°'|, of the leading jet for events with at least one jet with 'p'^ > 30 GeV and |j/J°*| < 4.4 
in the final state and (b) as a function of lyj"^*! of the second leading jet for events with at least two 
jets. The cross sections are normalized to the inclusive Z (— >■ it) cross section. The other details 
are as in Figure 2. 

increasingly from the data for larger values of K'^ and S't, which confirms and extends the 
results in reference [47] to a higher energy regime. The discrepancy is attributed to missing 
higher jet multiplicities in the fixed-order calculation. This interpretation is investigated 
further in what follows. 

Figure 16(a) shows, at reconstruction level, the average jet multiplicity as a function 
of H'x for the Z — )■ ee channel. Compatible results have been obtained in the muon chan- 
nel. Predictions by ALPGEN and SHERPA are consistent with the data. For values of 
i?T ~ 350 GeV, where data and NLO calculation start to deviate significantly, the average 
jet multiplicity exceeds two. A similar measurement is performed as a function of the p^ 
in the Z — ^ ///.t channel and shown in figure 16(b). Compatible results have been obtained 
in the electron channel. For values of p^ ~ 200 GeV, where the NLO predictions underes- 
timate the measured cross section (see figure 8), on average two jets are resolved, typically 
one hard jet that carries most of the Z recoil, accompanied by a soft jet. In both cases, the 
kinematic regions where the NLO fixed-order calculations perform poorly are characterized 
by average jet multiplicities in excess of the fixed order used in the NLO calculation. 

Figures 17(a) and 17(b) replace the fixed-order BlackHat+SHERPA estimate for H'x 
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Figure 10. (a) Measured cross section for Z (— >■ ££) + jets as a function of the absolute value of the 
rapidity, |y^°'|, of the third jet for events with at least three jets with pl^ > 30 GeV and |j/J°*| < 4.4 
in the final state and (b) as a function of \y'°^\ of the fourth jet for events with at least four jets. 
The cross sections are normalized to the inclusive Z (— >■ i£) cross section. The other details are as 
in Figure 2. 

and p^ in figures 15 and 8 with the 'exclusive sum' of the cross sections for the first two 
jets: (Z (+1 jet)) + (Z (+ > 2 jets)). The exclusive sum is consistent with the observed Ht 
and p^ spectra in the phase space considered. These results support the interpretation of 
the poor performance of the fixed-order calculation for inclusive quantities like H^, St and 
p^ as a sign of missing higher jet multiplicities. Agreement with the data can be restored 
by adding explicitly higher jet multiplicities via exclusive sums. 

11 Conclusions 

Cross sections for jets produced in association with a Z boson have been measured in 
proton-proton collisions at ^/s = 7 TeV with 4.6 fb^^ of data observed with the ATLAS 
detector at the LHC, using electron and muon decay modes of the Z boson. The data 
have been unfolded to the particle level and compared with predictions from the SHERPA 
generator, from MC@NLO interfaced with HERWIG, from the ALPGEN generator, inter- 
faced with HERWIG, and with fixed-order calculations from BlackHat+SHERPA. The 
cross sections are quoted with respect to a phase-space region defined by Z candidates 
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Figure 11. (a) Measured cross section for Z (— > tt) + jets as a function of the separation in rapidity, 
jAy-'^l, between the two leading jets and (b) as a function of the invariant mass of the two leading 
jets, m?^ ^ for events with at least two jets with j}!^ > 30 GeV and \y^'^^\ < 4: A in the final state. 
The cross sections are normalized to the inclusive Z (— >■ i£) cross section. The other details are as 
in Figure 2. 

constructed from opposite-sign leptons with j»t > 20 GeV, |?7| < 2.5, AR > 0.2 and 
66 GeV < m" < 116 GeV and for jets with pj^* > 30 GeV, |yj^*| < 4.4 and AR'^^ > 0.5. 

Cross sections as a function of the inclusive and exclusive jet multiplicities and their 
ratios have been compared, as well as differential cross sections as a function of transverse 
momenta and rapidity of the jets, angular separation between the leading jets and the 
inclusive variables ifx and 5t. Compared with previous publications, the sensitivity has 
been extended to regimes with larger jet multiplicities and larger jet transverse momenta. 
In addition, the sample has been compared to theory in specific kinematic regions governed 
by large logarithmic corrections. 

In general, the predictions of the matrix element plus parton shower generators and 
the fixed-order calculations are consistent with the measured values over a large kinematic 
range. MC@NLO fails to model not only higher jet multiplicities but also the transverse 
momentum of the leading jet. The transition from staircase to Poisson scaling of the exclu- 
sive jet multiplicity ratio, expected from theory when introducing a large scale difference, 
is observed in the data. 

In events where two jets have passed a VBF preselection, the cross sections for higher 
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Figure 12. (a) Measured cross section for Z (— > i€) + jets as a function of the distance in between 
the two leading jets, A(/)-'^ , and (b) as a function of the distance Ai?^^ between the two leading jets, 
for events with at least two jets with 'p'^ > 30 GeV and |yj'^'| < 4.4 in the final state. The cross 
sections are normalized to the inclusive Z (— >■ it) cross section. The other details are as in Figure 2. 



jet multiplicities are overestimated by ALPGEN. This leads to a small underestimation of 
the probability for Z+jets events to survive a veto on a soft third jet. 

ALPGEN predicts a too hard spectrum of the transverse momentum of the leading jet, 
of p^, i?T and S't in a regime where large corrections from higher-order electroweak and 
higher-order QCD processes are expected. The jet rapidity distribution is predicted to be 
too wide in BlackHat+SHERPA and in SHERPA. BlackHat+SHERPA underestimates 
the cross section for large p^ where more than one jet can be resolved. The H^ or 5t 
spectra predicted by BlackHat+SHERPA fixed order NLO calculations deviate by several 
standard deviations from the measured spectra in the hard H^ and S't regime characterized 
by large average jet multiplicities. The observed spectra of iJx and pf can be described by 
an exclusive sum of BlackHat+SHERPA fixed-order calculations for Z (+1 parton) and 
Z (+ > 2 partons). 
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Figure 13. (a) Measured cross section for Z (— > ti) +jets as a function of the transverse momentum, 
p^r^ , of the third jet and (b) as a function of the absolute value of the rapidity, |y^°'|, of the third 
jet, in events passing the VBF preselection (at least two jets with j}^ > 30 GeV and |2/J°*| < 4.4 
and m^^ > 350 GeV and [Ay^-'l > 3.0 for the two leading jets). The cross sections are normalized 
to the inclusive Z (— >■ £i) cross section. The other details are as in Figure 2. 
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Figure 14. Fraction of events that pass a veto on a central (|?7| < 2.4) third jet after VBF 
preselection (at least two jets with pi^ > 30 GeV and \y^'^*'\ < 4.4, m^-' > 350 GeV and 
\Ay^^ I > 3.0 for the two leading jets) as a function of the third jet pi^ threshold, min pl^ , (a) in the 
electron channel and (b) in the niuon channel, measured in data and predicted by the generators 
ALPGEN and SHERPA (see legend for details) . The data points indicate the measured distribution 
after subtraction of electroweak and multi-jet background. The hatched bands correspond to the 
combined statistical and systematic uncertainty on the Z + jets prediction, using ALPGEN to derive 
the systematic uncertainties. The error bars on each data point show the combined statistical and 
systematic uncertainty on the data. The bottom panel shows the MC/data ratio. The shaded band 
corresponds to the total systematic uncertainty and the error bars to the statistical uncertainty on 
the MC/data ratio. 
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Figure 15. (a) Measured cross section for Z (— )■ €^)+jets as a function of the scalar px sum of 
the leptons and the jets, i/x, and (b) as a function of the scalar px sum of the jets, S't, in events 
with at least one jet with 'p'^ > 30 GeV and ly-'"'! < 4.4 in the final state. The cross sections are 
normalized to the inclusive Z {—> ii) cross section. The other details are as in Figure 2. 
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Figure 16. (a) Average number of jets, <iVjet>, in Z (— >■ ee) +jets events as a function of the scalar 
Px sum of the leptons and the jets, iJx, and (b) average number of jets in Z (— >■ /x/u) + jets events as a 
function of the transverse momentum of the Z boson candidate, p^, measured in data and predicted 
by the generators ALPGEN and SHERPA (see legend for details). The data points indicate the 
measured distribution after subtraction of electroweak and multi-jet background. The hatched 
band corresponds to the combined statistical and systematic uncertainty on the Z+jets prediction, 
modelled with ALPGEN. The error bars on each data point show the combined statistical and 
systematic uncertainty on the data. The bottom panel shows the MC/data ratio. The shaded band 
corresponds to the total systematic uncertainty and the error bars to the statistical uncertainty on 
the MC/data ratio. 
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Figure 17. (a) Measured cross section for Z {-^ ££) + jets as a function of the scalar px sum of the 
leptons and the jets, H^, and (b) as a function of the transverse momentum of the Z candidate, p^, 
in events with at least one jet with pf^ > 30 GeV and jyj''* | < 4.4 in the final state. The cross sections 
are normalized to the inclusive Z (— >■ i£) cross section. The unfolded data are compared to NLO 
pQCD predictions from BlackHat+SHERPA, obtained by adding the exclusive Z (^ ££) +1 jet 
and the inclusive Z (— >■ ££)+ > 2 jets calculations and corrected to the particle level. The error bars 
indicate the statistical uncertainty on the data, and the hatched (shaded) bands the statistical and 
systematic uncertainties on data (prediction) added in quadrature. 
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Incl. jet multiplicity 


Data cross-section (pb) 


^had 


^QED 


> 1 jets 

> 2 jets 

> 3 jets 
>4 jets 

> 5 jets 

> 6 jets 

> 7 jets 


[ 6.88± 0.01 (stat)± 0.52 (syst) ± 0.13 (lumi) ]xlO^ 
[ 1.51± 0.01 (stat)± 0.15 (syst) ± 0.03 (lumi) jxlO^ 

3.09± 0.03 (stat)± 0.40 (syst) ± 0.06 (lumi) 
[ 6.55± 0.13 (stat)± 1.06 (syst) ± 0.12 (lumi) ]xlO"^ 
[ 1.35± 0.06 (stat)± 0.27 (syst) ± 0.02 (lumi) ]xlO"^ 
[ 2.53± 0.27 (stat)± 0.60 (syst) ± 0.05 (lumi) ]xlO"^ 
[ 6.23± 1.46 (stat)± 2.14 (syst) ± 0.11 (lumi) ]xlO"^ 


1.027±0.015 
1.036±0.017 
1.031±0.033 
1.043±0.023 


0.976±0.005 
0.979±0.005 
0.980±0.005 
0.982±0.004 


Incl. jet multiplicity ratio 


Data cross-section ratio 


^had 


^QED 


> 1 jets / > jets 

> 2 jets / > 1 jets 

> 3 jets / > 2 jets 

> 4 jets / > 3 jets 

> 5 jets / > 4 jets 

> 6 jets / > 5 jets 

> 7 jets / > 6 jets 


[ 1.42± 0.00 (stat)± 0.11 (syst) ]xlO"^ 
[ 2.18± 0.01 (stat)± 0.07 (syst) ]xl0"^ 
[ 2.05± 0.02 (stat)± 0.07 (syst) ]xl0"^ 
[ 2.12± 0.04 (stat)± 0.08 (syst) ]xl0"^ 
[ 2.06± 0.08 (stat)± 0.10 (syst) ]xl0"^ 
[ 1.87± 0.18 (stat)± 0.13 (syst) ]xl0"^ 
[ 2.46± 0.49 (stat)± 0.39 (syst) ]xl0"^ 


1.036±0.015 
1.009±0.002 
0.995±0.016 
l.OlliO.OlO 


0.995±0.010 
1.003±0.010 
l.OOliO.OlO 
1.001±0.009 


Excl. jet multiplicity ratio 


Data cross-section ratio 


^had 


^QED 


1 jet / jet 

2 jets / 1 jet 

3 jets / 2 jets 

4 jets / 3 jets 

5 jets / 4 jets 

6 jets / 5 jets 

7 jets / 6 jets 


[ 1.29± 0.00 (stat)± 0.10 (syst) ]xl0"^ 
[ 2.23± 0.01 (stat)± 0.08 (syst) ]xl0"^ 
[ 2.03± 0.02 (stat)± 0.07 (syst) ]xl0"^ 
[ 2.14± 0.05 (stat)± 0.08 (syst) ]xl0"^ 
[ 2.11± 0.11 (stat)± 0.10 (syst) ]xl0"^ 
[ 1.74± 0.22 (stat)± 0.10 (syst) ]xl0"^ 
[ 2.60± 0.79 (stat)± 0.45 (syst) ]xl0"^ 


1.032±0.013 
1.013±0.010 
0.990±0.032 
1.022±0.028 


0.994±0.010 
1.003±0.010 
l.OOliO.OlO 
1.001±0.009 


Excl. jet multiplicity ratio 
pj^* (1st jet) > 150 GeV 


Data cross-section ratio 


^had 


^QED 


2 jets / 1 jet 

3 jets / 2 jets 

4 jets / 3 jets 

5 jets / 4 jets 

6 jets / 5 jets 

7 jets / 6 jets 


1.04± 0.03 (stat)± 0.03 (syst) 
[ 4.82± 0.13 (stat)± 0.16 (syst) ]xl0"^ 
[ 3.71± 0.17 (stat)± 0.16 (syst) ]xl0"^ 
[ 2.85± 0.21 (stat)± 0.12 (syst) ]xl0"^ 
[ 2.67± 0.37 (stat)± 0.21 (syst) ]xl0"^ 
[ 2.57± 0.78 (stat)± 0.51 (syst) ]xl0"^ 


1.004±0.002 
0.989±0.037 
1.025±0.040 


1.000±0.009 
1.002±0.006 
0.996±0.006 



Table 4. Combined inclusive Z {-^ ^^)+jets cross sections per lepton flavour and their ratios 
and exclusive cross-section ratios for various preselections measured in data together with the 
corresponding non-perturbative corrections S^^"^ and (5^^^. The cross sections are quoted with 
respect to a phase-space region defined by Z candidates constructed from opposite-sign leptons 
with pT > 20 GeV, |?7| < 2.5, AR" > 0.2 and 66 GeV < m" < 116 GeV and for jets with 
pi;?* > 30 GeV, lyj'^*! < 4.4 and AR'^i > 0.5. 
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